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Hindered magnetic dipole transitions Y(nS) — > jr/b(n'S) are studied in the covariant light- 
front approach. Compared with the allowed magnetic dipole transitions, we find that results 
for hindered magnetic dipole transitions are sensitive to heavy quark mass and shape pa- 
rameters of the light-front wave functions. It is possible to tune the parameters so that the 
predictions of branching fractions of T(2S, 3S) — > 777;, are consistent with the recent experi- 
mental data, but the relevant decay constant of r\b is much smaller than that of T(1S). We 
also generalize the investigation to the charmonium sector and find the the same conclusion. 

<N 

QJ ' Since the photon energy in allowed magnetic dipole (Ml) transitions is limited, hindered Ml 

transitions with changes of principal quantum numbers are believed to offer better chances to 
discover the pseudoscalar quarkoniums such as 7/f,. Several years ago, the CLEO collaboration set 
the upper bounds for this kind of decays [l| 

B(T(35) -> 7%) < 4.3 x 1(T 4 , 

B(T(2S) -> 7%) < 5.1 x 10" 4 . (1) 



0,; 



In 2008 the BaBar collaboration observed a peak in the photon energy spectrum at E~ = 
(921.2^2 8 i 2-4) MeV in radiative Y(3»S) decay 0|. This is viewed as the first observation of the ryt, 
■ meson which is the lowest bottomonium ground state. The branching ratio of this radiative decay 

G\ '. i s 

t> 
in 

m 1 B(T(3S) -> 7r/ 6 ) = (4.8 ± 0.5 ± 1.2) x 10~ 4 . (2) 

tN 

CO The measured mass 

o ' 

m(jfo) = (9388.9i|^ ± 2.7)MeV 
^ I gives the mass split 

5m %(ls) = m(T(15)) - m( Vb ) = (71.4±|? ± 2.7)MeV. (3) 
Subsequently the % meson has also been observed in the radiative T(25) decay 0] 

B(T(2S) -> 7%) = (3.9 ± l.ltH) x lO" 4 . (4) 
The ratio of the two branching fractions is measured as 

D _£(T(2S)^ 7 %) no _ 

K = B(T(35) 7%) ~ a24 -o.i9- (5) 



The recently updated results by the CLEO collaboration [4| are consistent with measurements by 
the BaBar collaboration taking into account the uncertainties 

B(T(3S) -> 7%) = (7.1 ± 1.8 ± 1.1) x 10" 4 , (6) 
B(T(25) -> 7%) < 8.4 x 10" 4 . (7) 
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Decay widths of Ml transitions can be expressed in a well-known formula 

2 

(8) 



T(n 3 5i ->■ 7n /1 5o) = ^aen-^- 
S ^ mi 
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a = 1/137 is the fine-structure constant, 6q denotes the charge in unit of |e| of the transition 
quark, and thq denotes the quark mass, k is the energy of the photon in the vector meson rest 
frame. For allowed transitions (n = n'), the emitted photon is typically soft. One can expand the 
spherical Bessel function j'o(if) by k. Radiative corrections and nonrelativistic corrections can be 
systematically studied in the effective field theory 0]. In particular the predicted branching ratio 
of J /if) —7- 777c is well consistent with the experimental data. 

The emitted photon in hindered Ml transitions is rather energetic and the convergence of the 
expansion in k becomes poor. On the other aspect, the decay amplitude is zero at the leading 
power in k since the different wave functions are orthogonal. Higher power terms contribute and 
they are sensitive to the treatments of the nonperturbative dynamics in transition form factors. For 
more than twenty years, hindered Ml transitions have received extensive interests but theoretical 
predictions vary over several orders 0| . The recent experimental data could offer an opportunity to 
investigate the dynamics in the quarkonia transition and in particular might be helpful to constrain 
the form of the wave functions. 

In this work, we use the covariant light-front quark model (LFQM) to investigate the Ml 
transitions and examine whether the commonly-used light-front wave functions (LFWF) can si- 
multaneously explain the experimental data for both allowed and hindered Ml transition. The 
light-front QCD may be the only potential candidate to reconcile the low energy quark model and 
the high energy parton model. The LFQM Qj can give a full treatment on spins of hadrons using 
the so-called Melosh transformation. Physical quantities are represented as the overlap of LFWF. 
These wave functions are expressed in terms of the internal variables of the quark and gluon degrees 
and thus are manifestly Lorentz invariant. To preserve the covariance and remove the dependence 
of physical quantities on the direction of the light-front, Jaus proposed the covariant LFQM in 
which the zero-mode contributions are systematically included [8(. The application of the covari- 
ant LFQM to the s-wave and p-wave decay constants and form factors is very successful [ SMl4l|. 
Under this framework, the transition T(15) — > jr]b(lS) has already been studied in Ref. [la, [la] 
(see also Ref. 0). 

In the following we employ the light-front decomposition of the momentum P = (P~ , P + , P±), 
where P ± = P° ±P 3 . These momenta can be expressed in terms of the internal variables (xi,p±) 
as: 

Pl,2 = %l,2P + ,Pl,2± = x 1)2 P± ± p±, (9) 
with x± + x 2 = 1. With these internal variables, one can define some useful quantities 



X\ X2 



_ X 2 M m|+pj _ / 2 + „2 +V 2 (w) 

here ej can be interpreted as the energy of the quark/antiquark and Mq can be viewed as kinetic 
invariant mass of the meson system. The transition amplitude of V(P, ey) — > P / y*(q, e 7 ) is usually 
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parametrized as 

A = iee^ pa eye v q p P°V{q 2 ), (H) 

where the photon is firstly taken as off-shell q 2 7^ and the convention e 0123 = 1 is adopted. 
Feynman diagrams for this process are given in Fig. [H in which the photon is emitted from the 
quark shown in (a) or from the antiquark shown in (b). Due to the charge conjugation invariance, 
these two diagrams provide identical contributions. It is straightforward to evaluate these two 
diagrams: 

A - iee Q l\ c J \ NiNzN ^ V + NxN2 N' 2 s ^} e 7 e v> W 



where iVj = pf — m\, N- = p' 2 — m 2 and 

Plu ~ P2u 



s a = Tr 



= Tr 



Wy 



{-fa + m 2 )^{A + m 1 )7 M (^i + mi) 



y u ~ ^Wv^) ^ 2 + m2)1 ^-A + ^2)75(^1 + rnx) 



(13) 
(14) 



Functions Hy, Hp, depending on the four momentum of the internal quarks, are the wave functions 
for vector and pseudoscalar mesons. In the absence of singularity in the Hy, Hp, integrating over 
the minus component of the internal momentum will pick up the pole in the propagators. Then the 
functions Hy p are replaced by hy,p which only depends on the plus component and the transverse 
momentum and Wy is replaced by wy 



wy = M + mi + rri2, (15) 

where (t>(x2,p±) denotes the momentum distribution inside the meson. In order to preserve the 
covariance, it is suggested that the zero mode contributions should be added in the covariant 
LFQM. The inclusion of the zero-mode contribution corresponds to a proper way to integrate out 
the minus component and after the integration expressions for the form factor V{q 2 ) are given as 



follows [la. lid] 



v{q) = 8^J dX2d p± — ^mok — r Q + + ~~^~ ) ) ■ (16) 

and the decay width of V — > P7 is then evaluated as 

F(V -> P 7 ) = l^ K-^) 3 |y(Q)|2 , (17) 

In the niQ — > 00 limit, the heavy quark inside the quarkonium moves nonrelativistically. In- 
tegrating out the hard off-shell degrees of freedom, one reaches an effective field theory known as 



nonrelativistic QCD [18(] to deal with the low energy dynamics. In a heavy quarkonium, the square 
of the transverse momentum p\ is of the order Aq CD and the momentum fraction of the heavy 
quark x 2 — 1/2. The kinetic invariant mass {Mq,M'q) is roughly 2mQ. At the leading power in 
1/rriQ, the form factor V^(0) in Eq. (fl~6j) is reduced to 

V (°) - l dx 2 d Pi- — ■ ( 18 ) 

87^ J m Q 
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FIG. 1: Feynman diagrams for V — > P7. 



This formula, with the LFWF discussed in the following, could formally reproduce the leading 
power behavior shown in Eq. ([8]) in the momentum space. 

For hindered Ml transition, the transition amplitude is from the higher power corrections and 
one needs a specific approach to calculate them. As we can see from the above discussion, the 
form factors in the covariant LFQM are expressed as the overlap of the LFWF. It contains the 
nonperturbaitve coupling of the quark pair with a meson. Except for some limited cases, these 
quantities can not be derived from the first principle. Practically one often adopts a phenomenolog- 
ical form. The commonly-used wave function for the 5-wave mesons with lowest principle number 
is the Gaussian-type 



<t>(%2, P±) 

dp z 
dx 2 



P 2 
eie 2 



3/4 



dp z 
dx2 



exp 



pj +Pl 
2/3 2 



(19) 



xix 2 M 

In the case of the wave functions of the 25, 3S mesons, we adopt the solutions of the harmonic 
oscillator problem 

3/4 /3 

U~~J exp 

3/4 

exp 



<hs{r) 
<t>zs(r) 



t)' 



/3V 







/3V 
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^ - 5/3 V + /3V 



(20) 



The Fourier transformations, with the Jacobi determinant V^n 2 -, are derived as 



4>2S(X2,P±) 
4>3S{X2,P±) 



VI 
VI 



3/4 



dp z 
dx 2 



exp 



3/4 



dp z 
dx 2 



exp 



P 2 ±+P 2 z 

2/5 2 

P±+P 2 z 

2/3 2 



P±+Pj _ 

{pI+pI? 

(3* 



HpI+pD 
p 2 



This kind of LFWF for the 2S mesons has been used in Ref. Uj while its light-cone form is studied 
in Ref. 121]. It would be meaningful to use the allowed and hindered Ml transitions to constrain the 
forms of the LFWF and in particular to test the harmonic-oscillator wave functions given above, 
as most of the previous studies in the covariant LFQM mainly focus on mesons without any radial 
excitation. 

Two kinds of inputs are required in the numerical analysis: masses (for hadrons and constituent 
quarks) and the shape parameters in the wave functions. For the heavy quarks, the difference 
between the current mass and the constituent mass is small. From the PDG 20(], the current 
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TABLE I: Input parameters in the light- front wave functions (in units of GeV). 





T(15) 


T(25) T(35) 


A /Til 

Mass 


9.46030 


10.02326 10.3552 


/t 


0.7150 ± 0.0048 0.4974 ± 0.0045 0.4301 ± 0.0039 




1 QOQ+0.037 
l.OZO_ q 35 


o-94D±g:g2 o.8oai8:S 




JH> 


V>(2S) 


Mass 


3.069616 


3.68609 


fi> 


0.4163 ± 0.0053 0.2961 ± 0.0025 





641 + 027 
u.o^i_ 025 


n 4Q7+0 019 

U -^ J ' -0.017 



masses are m r = (1.27toHl) GeV and m b = (4.20^;^) GeV in the MS renormalization scheme. In 
this work, we will choose three different values for the constituent quark masses 

m c = (1.1, 1.3, 1.5) GeV, m b = (4.0, 4.4, 4.8)GeV. (22) 

As for hadron masses, all of them used in the present analysis have been measured except the one 



of r/b(2S'). We will use the mass split in the charm sector [20(] 



Sm Vc(ls) = 116.5MeV, 5m Vc(2S) = 48.1MeV 



to estimate the mass split for the r]b(2S) 

<K»(2S) = s mVb(1S) Sm Vc(2S) = 29.5MeV. (23) 

Shape parameters of the LFWF are usually determined by decay constants of hadrons whose 
expressions are given in Ref. 0]. Decay constants of T(nS) are extracted from the partial decay 
width of the leptonic T — > e + e~ decays 

r. = r(T-» e + e -)-*!^. (24) 

The shape parameters (corresponding to m;, = 4.4 ±0.4 GeV and m c = 1.3 ±0.2 GeV) are collected 
in tablelH In the case of % mesons, their decay constants are not measured at present. In the future, 
these decays may not be well constrained as the uncertainties in the ideal mode r)b — > 27 would be 
typically large. In this work, we first try to tune the shape parameter so that the predictions of 
branching fractions of T(25', 35) — > 7% decays are consistent with the data provided by the BaBar 
collaboration. 

With above parameters, results for the branching ratios in T decays are collected in table [II] 
and table IIIIl in which we have used three sets inputs for the shape parameters and three different 
values for the quark mass. Several remarks are give in order. Firstly, from these tables one can 
see that the decay constant of the r] b is not very sensitive to the shape parameter and the quark 
mass m b - Secondly the allowed channels T(nS) — > rj b {nS)(n = 1,2) are not sensitive to these 
inputs either. On the contrary, results for hindered channels are strongly dependent on the input 
parameters, which is not beyond expectation. If we adopt the mj = 4.4 GeV and f3 Vb = 0.93 
GeV, theoretical predictions of branching fractions of T(25, 35) — > 7% are consistent with the 
experimental data given in Ref. 0,0]. The ratio of branching fractions is predicted as 0.66. 



TABLE II: Branching ratios (in units of 10 4 ) of Ml transitions: Y(nS) — > 777b (n'S). The b quark mass is 
used as mj, = 4.4 GeV. 







0.90 0.93 


0.96 






389 406 


424 


B(T(15) 


->%) 


2.7 2.8 


2.9 


B(T(25) 


-> %) 


21.3 5.0 


0.003 


B(T(3S) 


->»».) 


2.3 7.7 


18.6 


R 




9.4 0.66 


1.5 x 10~ 4 




0.75 0.8 


0.85 




343 372 


400 


B(T(2S) -» 


»7b(25)) 


0.3 0.3 


0.4 


B(T(3S) -»■ 


7?6(2S')) 40.3 3.1 


6.7 



TABLE III: Branching ratios (in units of 10 4 ) of Ml transitions: Y(nS) —> 7776 (n'S). The shape parameters 
of rjb and rjb(2S) are used as (3 Vb = 0.93 GeV; P nb (2S) = 0-80 GeV. 







4.0 


4.4 


4.8 






419 


406 


394 


B(T(1S) 


%) 


3.3 


2.8 


2.5 


B(T(2S) 


77b) 


24.8 


5.0 


0.003 


B(T(3S) 


77b) 


3.3 


7.7 


14.4 


R 




7.5 


0.66 0.0002 


m b 




4.0 


4.4 


4.8 


fvb(2S) 


379 


372 


364 


B(T(25) -> 


Vb (25)) 


0.4 


0.3 


0.3 


S(T(3S) -> 


77 b (2S)) 18.2 


3.1 


0.055 



In the heavy quark limit, pseudoscalar mesons and vector mesons belong to the same super- 
multiplet and their decay constants are almost the same. If we adopt the same decay constant for 
rjb with its vector partner, the shape parameter is determined as /5L = 1-451 GeV and theoretical 
predictions of branching ratios are given as 

B(T(15) -> 7776) = 3.2 x 1(T 4 , 
B(T(2S) -> 7%) = 5.8%, 

B(T(35) -)■ 77/fe) = 13.7%. (25) 

The latter two results are much larger than the experimental data given earlier. These inconsis- 
tencies between theoretical results and the data may imply either the inappropriate form of the 
LFWF used in this analysis, or a small decay constant for r/ b compared with that for T. 

As another check, we will extend the above discussion into the charmonium sector: ip(nS) — > 



"fr] c (n'S). The relevant experimental data is given as 20. 1211] 



B(J/tP -> 77?c) = (1.7 ±0.4)%, (26) 
B(ip(2S) -»• 777 C ) = (3.4 ± 0.5) x 10" 3 , (27) 
B(ip(2S) -»• 77? C (2S)) < 7.4 x 10" 4 . (28) 
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TABLE IV: Results for the branching ratios of Ml transitions involving r\ c in the covariant LFQM (in units 
of 10~ 2 for J/ijj — > r] c and 10 -3 for ip(2S) — > r) c ). In the left sector, the mass is used as m c = 1.3 GeV; while 
in the right sector, the shape parameter is /3 Vc = 0.55 GeV. 



Vc 0.5 0.55 0.6 


m c 1.1 1.3 1.5 


f Vc 253 281 308 


f Ve 283 281 276 


B{J/ip^r lc ) 3.0 3.1 3.1 


B(J/i/)^r) c ) 3.8 3.1 2.5 


B(ip(2S) -> r, c ) 5.5 4.2 33.4 


B(tp(2S) -> i7 c ) 0.6 4.2 15.3 



Again with proper inputs the results in table IIVI could be consistent with these data, but the 
relevant decay constant of rj c is much smaller that of its vector partner J/tp. Our result for 
Bty(2S) -> r, c (2S)) 



BmS) -> 1Vc (2S)) = (6.33:£i$ x 10" 4 . (29) 

is approaching the experimental upper bound. The first uncertainties are from m c = (1.3±0.2)GeV 
and the second uncertainties are from fi^jiS) = (0.45±0.05) GeV. In Ref. [16(, the inputs m c = 1.56 



GeV, Pr) c = 0.820 GeV and /3j/^ = 0.613 GeV are adopted and prediction of B(J/ip(lS) — > jr] c ) is 



consistent with the data reported by the CLEO collaboration [221 ] . If we adopt these parameters, 



the branching ratio of ip(2S) — > jr] c is larger than the experimental data by two orders 

B{^{2S) -> 7r? c ) = 18.2%. (30) 

Inspired by these inconsistencies, one may get the same conclusion as in the bottomonium sector. 

In conclusion, we studied the hindered Ml transitions ip(2S) — > ji] c and the T(nS) —> 
777b 2S) in the covariant light-front quark model. Compared with the allowed Ml transitions, 
theoretical results for hindered channels are found to strongly depend on quark masses and the 
shape of the light-front wave functions. Because of this sensitivity, these transitions offer good 
laboratories to study the dynamics between the transitions and constrain the forms of the light- 
front wave functions. Using the harmonic-oscillator wave functions and the experimental data, our 
results show that the decay constants of the pseudoscalar mesons are required to be much smaller 
than the decay constants of their vector partners. With the precise decay constants measured in 
various other channels in the future, different forms of the wave functions could be directly tested 
in the channels studied here. A more practical way in the future is to incorporate more radiative 



transitions between heavy quarkonium states as in Ref. [23 1. 

Note added: During the long-time revision of this work, a similar paper [24[ | appears. One of 
the important differences is: we have shown that hindered Ml transitions are sensitive to the shape 
of the wave functions as they should. 
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